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Abstract: We study a meson mass splitting due to isospin violation in holographic dense 
matter. We work in a D4/D6/D6 model with two quark flavor branes to consider asym- 
metric dense matter in holographic QCD. We mainly consider two cases. We first consider 
rri^ jmT ~ m^jmu to study the effect of isospin violation on the meson masses. Then, 
we take rri^ jmr ~ mg/mq, where niq ~ [rriu + m(i)/2, to calculate in-medium kaon-like 
meson masses. In both cases we observe that the mass splitting of charged mesons occurs 
at low densities due to the asymmetry, while at high densities their masses become de- 
generate. At intermediate densities, we flnd an exotic behavior in masses which could be 
partly understood in a simple picture based on the Pauli exclusion principle. 
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1. Introduction 

Asymmetric dense matter is a fascinating subject in modern nuclear physics. Nature 
provides isospin asymmetric nuclear matter in the interior of a nucleus with unequal number 
of protons and neutrons and in compact stars such as a neutron star. Through heavy ion 
collisions with neutron rich, stable and/or radioactive nuclei, we could study the properties 
of this matter in the laboratory. As density goes up, strangeness starts to play important 
roles. For instance, the transition from nuclear matter to strange matter, through onset of 
hyperon matter or kaon condensation, shows a critical impact on the properties of compact 
stars such as the mass-radius relation and on core-collapse supernovae. 

In addition to vacuum properties of dense matter, the properties of hadrons in such 
environment are of interest. The in-medium mass of K~ is one of the key quantities to 
determine the onset of kaon condensation. The in- medium modification of kaon/antikaon 
properties such as the mass can be observed experimentally. One of the experimental 
observables is the K^/K^ production ratios. Upcoming experiments such as the GSI 
Facility for Antiproton and Ion Research (FAIR) would reveal the in-medium properties of 
heavy-light mesons, D and D, through experimental signatures of their production ratio 
and J/ip suppression. We refer to [|| for a review on the issues in dense matter. 

The main goal of this paper is to study the in-medium meson masses and their mass 
splitting in an asymmetric dense matter using a holographic QCD model. Based on the 
AdS/CFT correspondence 0, a holographic model of QCD [Q, see [11] for a review, is 



offering a new analytic tool to study dense matter. In g| , a simple model for nuclear matter 
to strange matter transition has been proposed in a D4/D6/D6 model. In this work the 
authors introduced two fiavor D6 branes for light (u or d) and intermediate (strange) quarks 
respectively and compact D4 brane with Nc fundamental strings attached. By considering 



- 1 - 



the force-balancing condition and energy minimization, they could study transition from a 
dense matter with only light quarks to a matter with both light and intermediate quarks. 
Since the background in is confining D4 geometry, the authors could associate this 
transition with nuclear to strange matter transition. Our present study is based on this 
D4/D6/D6 model to gain asymmetry in dense nuclear matter. In this sense, the mass 
difference of two quarks is the origin of our asymmetry. To be more realistic, our model 
should be supplied with a few more ingredients. For instance, in neutron stars the charge 
neutrality favors neutrons rather than protons since protons should come in together with 
electrons. A real neutron star contains a small fraction of protons and electrons to prevent 
neutrons from undergoing a weak decay into protons and electrons, though. Obviously we 
need to have three flavors with different quark masses. To this end we have to extend the 
present model to a D4/D6/D6/D6 model. With these cautions in mind, we will study the 
meson masses in asymmetric dense matter, hoping that our study would catch some main 
features of the physics of asymmetric matter. We will consider two cases with different 
quark mass ratios: m~^/m~ ~ rUii/mu and rn^ jmT ~ ms/mq, where ruq = {niu + m(i)/2. 
The first case is for a nuclear matter with isospin violation, and the second one is to mimic 
a transition from nuclear matter to strange (hyperon) matter. 

In Appendix we use a different value of the 't Hooft coupling from the one used in 
the main text to see if our results are sensitive to the choice of the coupling. 

2. Asymmetric dense matter 

In this section we review the work of Q, where a toy model, based on a D4/D6/D6 setup, 
for nuclear matter to strange matter transition has been proposed. In HI, the matter with 
one light and one intermediate quarks is named as strange matter without proper discussion 
of the strangeness quantum number. This is basic reason why the authors of Q considered 
their model as a toy. Note, however, that in some cases describing strangeness physics 
with two flavors works well. For instance, the essential physics of kaon condensation could 
be captured in the V-spin approach, where chiral SU(3)xSU(3) is effectively reduced to 
SU(2)x SU(2) ||2^. In this work we will make use of the asymmetry obtained in Q to 
study meson mass in asymmetric dense matter. 

Now we review the paper |[. The non-supersymmetric geometry for confining 
background of D4 is given by 

ds' = (j^y {r^,.dxV + fiU)dxl) + (-^ + U'dnl^ 

e' = 9s{^y\ FA = '-^e,, fiU) = l-{^y, R^ = .g.NJl (2.1) 

The Kaluza-Klein mass is defined as inverse radius of the X4 direction: Mkk = § j^/i ■ 
The bulk parameters, UkKjOs, and R, and the parameters of the gauge theory, Mkk and 
9ym^ are related by 

<"= 2,lJ,M,, ' V.K = IXM,,P., IC = ^, X=,ImN.. (2.2) 
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To make transverse space to D4 brane be flat (up to overall factor), we introduce new 



coordinate ^, 



TPW) 



, and obtain, in Euclidean signature, 



rr\ 3/2 



dx^ + /(C/)dxl) + 



3/2 



u 



(dp + f dfll) 



(2.3) 



The relation between U and ^ is 



jj3/2 _ ^/2 -j^ _|_ 



(2.4) 



where U' 



3/2 



^^"^ rescale ^ coordinate such that singularity is located at ^ = 1. 



With these new coordinates we rewrite the background ( |2.3p and the dilaton as 

/ t \ 3/2 1/3 

ds' = J {dt^ + + /(^)dxi) + (i?3e^^) ^ (dp + e^dQi) , 

e<^=(^l^y^'p/4^V2. (2.5) 

A baryon in four-dimensional theory could be described by a D4 brane wrapping S^. 
In this configuration, the background four-form field strength couples to the world-volume 
gauge field A(i-^ via the Chern-Simons interaction. The source of these world-volume gauge 
field can be interpreted as the end point of fundamental strings. 



We rewrite again the background metric (2^) as 

3/2 



ds'=(^j^^ e'^u+{de+dx^ + 



f{U)dx: 



+ {R'iKK) u:fe'^ [^+de^ + sin2 ed^i^ . (2.6) 

We take (t, 6*, 6a) as a world-volume coordinate of D4 brane, where Oa are angular coor- 
dinates on and turn on the U{1) gauge field on it. For simplicity, we assume that the 
position of D4 brane and the gauge field depends only on 9, i.e., ^ = S,{9), At = At{6), 
where 9 is the polar angle. The induced metric on the compact D4 brane is 



l + ^\d9'^ + sin^ 9dnl 



(2.7) 



Then the DBI action for the single D4 brane reads 

Sda = -//4 J e-'^Vdet(5 + 27ra'F) + ^4 j A G(4) 
= r4 y dtd9 sin^ 9 

dtCDA, 



(2. 
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where 



T4 = fJ.4^39s ^R^CkK, At 



2'Ka'At 



KK 



(2.9) 



Substituting the solution of equation of motion for gauge field At, we get 'Hamiltonian' of 
D4 brane as 



OF 



where 



uj do^Lo^fie + i'^)\lD{eY + sin6 e, 



D{e) = -2 + 3(cos e-\ cos^ e). 



(2.10) 



(2.11) 



Here, we set all fundamental strings are attached on north pole of D4 brane, because we 
are interested in the system in which fundamental strings connect baryon D4 brane and 
probe D6 branes as discussed in j 



After imposing boundary condition ^'(0) = and ^(0) = at = 0, we can get numer- 
ical solutions which are parameterized by initial value of ^o- The solutions corresponding 
different ^o's are drawn in Fig. 1(a). 





(a) 



(b) 



Figure 1: (a) Shape of D4 brane in 0) plane for different = 1.001, 1.05, 1.1, 1.2, 1.5 from 
inside to outside. Here south pole corresponds to 9 = {). (b) Force at the cusp of D4 brane for 
different ^o- 



If we denote the position of the cusp of D4 brane by Uc, the force at the cusp due to 
the D4 brane tension can be obtained by varying the Hamiltonian of D4 brane with respect 
to Uc while keeping other variables fixed; 



Fda 



dn 



dU, 



fix other values 
3^ 



NcTf 



1 + t 



72' 



(2.12) 
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where Tp = ^ is the tension of the fundamental string. As shown in Fig. 1(b), the 
tension at the cusp of D4 brane is always smaller than the tension of the Nc fundamental 
strings. Therefore, if there are no other objects, the cusp should be pulled up to infinity 
and the final configuration of D4 brane would be 'tube-like' shape as in |13]. To have a 
stable configuration, we put probe D6 branes which the other end point of the fundamental 
string can be attached. 

Now we consider the system with two fiavors, one light and one intermediate mass 
quarks, by introducing two D6 branes. To study dynamics of multi D-branes, we start 
from the non-Abelian DBI action 01 . 



5 = -/i6 J (fa STr[e-y-det(P[G,, + Gra{Q-^ - 6)''^Gbs] + T-^Frs)^detQ%\ , (2.13) 

where STr denotes the symmetrized trace for fiavor indices and P[ ] does the pull-back of 
lOD tensor to the world-volume of the branes. The matrix Q'*^ is defined as 

Q% = 6'l + iT[X'',X^]Gcb, (2.14) 

where T=l/{27ra') and X"'s are the coordinates being transverse to the branes which are 
2x2 matrix valued functions in general. 

To obtain the embedding solution of two D6 branes, we rewrite the bulk metric 

as 

3/2 



4; 

1/3 

+ {R^kk) ^ + p^dnl + dy^ + y^dcp^) , (2.15) 

where D6 brane world- volume coordinates are {t,x, p,9a)- We assume that the transverse 
coordinates depend only on p. The induced metric on D6 brane is 

ds^ = (k^Y%y^uj+ (dt' + dx^) 



R 



+ {R'iKKf" ^ [(1 + + y'<t>") dp" + p^dfii] , (2.16) 

where prime denotes to derivative with respect to p. Here, y and (p which transverse to D6 
brane world- volume are matrix- valued functions. 

In this work, we take the diagonal ansatz for embedding functions, which is physically 
relevant since we expect off-diagonal quark condensates to be zero, < ud >=< du >= 0. 

where is the position of each probe brane. On each brane, we turn on ^7(1) gauge field 
whose source is the end point of fundamental strings 
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In terms of the Pauli matrix 



10 1 
2 1 



1 /O - 

2 ii J 



1 n 1 

2 lo -Ij 



y can be written as 



where 



y = y^T^ + 



= y^ + y , y^ = y^ -y~ 



(2.19) 

(2.20) 
(2.21) 



Near the boundary, y^'^ becomes 



y'- 



m + m ) + 



y ~ 1 m — m 



+ c" 

c+ — c 



+ 



+ 



(2.22) 



If we put two branes on top of each other, the value of y^ becomes zero at the boundary. 
This is physically obvious since the exact two flavor isospin symmetry U(2) guarantees 
771+ = rn~ and c+ = c~. 

Similarly, U{1) gauge field can be written as 



At = A^T^ + A^T^. 

Near boundary, we can define isospin chemical potential and isospin charge as 



(2.23) 



A0 = (/.+ + ^-) + 
A' = (/x+-^-) + 



Q+ + Q- 
Q+ -Q- 



+ 



+ 



Qt 

fJ-T-\ n- + 

Pl + ^ + ■ 
P^ 



(2.24) 



With this diagonal ansatz, the commutators in the non-Abelian DBI action ( 2.13| ) 
vanish, and the action can be expressed as the sum of two independent embeddings 

Sd^ = -re j dtdpp^u:f^Ju:f{l + y'^)-F^ 

-TeJ dtdpp^ujf^J/\l + y'2) _ F2 ^ (2.25) 



where 



27ra'Ftp 



(2.26) 



Here it denotes that each DBI action is function of y^ and Af^ separately. We define a 
dimensionless quantity from the equation of motion for F^; 

p'^fF 



dSm 
dF^ 



Q 



± 



(2.27) 



± 
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It is related to the number of point sources (number of fundamental strings) by 



27ra'T6 

Through the Legendre transformation, we obtain the Hamiltonian 

T~iD& = F—S^ — Sd6 
oF 



(2.28) 



re / Wa.+ (y+)4/3 Q+2 + ^4^^(y+)8/3K/i + y+ 



.12 



+r6 / dpW^+(y-)4/3 Q-2+^4^^(y-)8/3U/l + y 



J2 



(2.29) 



Since the Hamiltonian of two D6 branes is split into two independent Hamiltonian of each 
brane, we can get the embedding solution of each brane separately. In our configuration, 
fundamental strings are attached on the upper brane and Q~ strings end on the lower 
brane, while the other end points of the fundamental strings are attached on the D4 brane. 
Therefore, we expect that two probe branes and spherical D4 brane meet at the same point 
eventually. The force at the cusp of the D6 branes is given by 



Fd6 



+ 



dUr. 



2na' V 1 - 6 



+ 



1 



Vc 



_/2 



27ra' V 1 - 6 



—3 



l + Vc 



_/2 



(2.30) 



^F+^{Q+) + F-^{Q~), 

where y^' denotes the slope of each brane at the cusp. To make the whole system stable, 
the force at the cusp of D4 brane should be balanced to those of D6 branes; 

Q 



Nr 



Fda = FUQ+) + F^JQ 



(2.31) 



where Q = + Q . Rewriting and yf by using new parameters a and /3, 

g+ = aQ , Q~ = (1 - a)Q, 



the force balance condition ( |2.31 ) becomes 



Vc, 



12 



+ 



(2.32) 



(2.33) 



Together with the force balancing condition, we consider the energy minimization of the 
system. The total energy of the system is 



(2.34) 
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where Ei is numerical integration of each 'Hamiltonian' in ( p.lCl| ) and ( 2.2E ) without overall 
constant T4 and tq. As an example, we fix the asymptotic values of probe branes as m"*" = 5 
and m~ = 0.1. The a dependence of the total energy with different densities is drawn in 
Fig. ^(a) and the density dependence of a for stable and minimum-energy configuration is 
shown in Fig. |^(b). In low density, only the lower brane (m~ = 0.1) constitutes a physical 
configuration. As density increases, at certain density, baryon vertex which attached to 
the upper brane can come into the system. And, at high densities the system consists of 
baryons attached to upper and lower branes equally. 




(a) 




Figure 2: (a) a dependence of regularized free energy, (b) Density dependence of a. 



We can also calculate the chemical potential difference of two quarks, m"*" and m", and 
its canonical conjugate density by using ( p. 24 ). Numerical results of density dependence of 
these quantities are drawn in Fig. ^. At low density, all fundamental strings are attached 
on the lower brane, i.e. = 0. In this regime, the density difference is nothing but 
Qi = —Q~, and it decreases linearly up to the transition density. Above the transition 
point, the number of fundamental strings on the upper brane increases so that the difference 
between and Q~ decreases. Therefore, the density difference increases and it approaches 
zero at large density as in Fig. ^(a). 

The density dependence of the chemical potential difference is shown in Fig. ^(b). 
The behavior at low density is similar to the density difference. However, at the transition 
density, the chemical potential difference suddenly becomes zero and keeps this value for 
all density afterwards. This behavior is quite reasonable in the following sense. A simple 
reason for the transition from a = to a 7^ is that the chemical potential (or Fermi 
energy) of light quark ^~ becomes comparable to the mass of heavier one m"*", /i^ ^ m~^, 
near the transition density. While, just above the transition point, the chemical potential 
of the heavier quark will be close to its mass, /x^ ~ . Therefore we could expect 
/x"*" — /U~ ~ at and just above the transition density. As density increases further, the 
chemical potential difference will be stay close to zero since at high densities the quark 
mass difference could be neglected. 
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(a) 



(b) 



Figure 3: (a) Density dependence of tlie density difference, (b) Density dependence of the chemical 
potential difference. 



3. Meson mass in asymmetric matter 

We begin this section with a simple example of isospin matter to gain a rough picture on 
the behavior of meson masses in asymmetric matter. Isospin matter with finite isospin 
chemical potential and zero baryon chemical potential was proposed as a useful system 
to improve our understanding of cold dense QCD [Q], though this kind of matter hardly 
exists in nature. At small (negative) isospin chemical potential fij with zero isospin number 
density, the mass of the pion is given by 

m^± = m^o + q\fii\ , (3.1) 

where q is the isospin charge of the particle. Note that vr'^ has zero isospin charge and so 
the mass will be intact with a finite isospin chemical potential. Since vr" and vr"*" have the 
opposite isospin charge q = —1 and q = +1 respectively, their mass will be splitting with 
increasing \fii\, linearly in this case. At some large value of the chemical potential \^J-i\, 
pion condensation will set in, and the in-medium behavior of the pion mass becomes more 
intriguing. For more details on isospin matter we refer to 

Now we calculate meson masses in asymmetric matter based on Q. To this end, we 
consider non-Abelian DBI proposed in Q. We will basically consider two cases. We first 
consider /m~ ~ md/rriu to study the effect of isospin violation on the meson masses. 
Then we will take rri^ jm^ ~ ms/mq, where niq = {rriu + m(^)/2, to calculate in-medium 
kaon-like meson masses. 

We remark that we are focusing on a general tendency of the in-medium meson mass 
not on the exact numbers since our model is yet not that close to a real QCD. 

3.1 D6 brane fluctuations 



We start from the non-Abelian DBI action ( 2.13| ) for D6-branes in the confining D4 back- 



ground. Since we are considering one light and one intermediate mass quarks, we take 
Nf = 2. 



S = -^i&j d'a STV [e-y-det(P[G„ + Gra{Q-^ - 5Y^Gbs]+T-^Frs)^J detQ\\ (3.2) 



For later convenience, we define some notations on the flavor matrix. In our case with two 
flavors, these matrices can be represented by 2 x 2-matrices 



M = Mqt^ + MiT' 



I Mq + Ms Ml- iM2 



'M+ Mi2^ 
1M21 M_j 



(3.3) 



2 \ Ml + iM2 Mo- Ms I 
using Pauli matrices. For diagonal brane embeddings the above action is expanded as 

-/i6 [(fa STv [e-*./-det{G,, + DrX'^D.XbiGab - iT[X'^, Xd]GacGdb) + T-iF^} 



S 



X (l + ^[x^x'']G5e[x^x'^]G, 



where the covariant derivative is defined by 

DsX^ = dsX'' + i[As,X^] . 
The induced metric on the D6 brane takes the form 



1/2 u;+ 



1/3 ^ 



(3.4) 
(3.5) 

1 + y'^) dp^ + p^dnf\ , (3.6) 



where y' = dy/dp. Embedding solutions for y are summarized in Sec. |2|. 

Now we consider some fluctuations to study meson spectrum in asymmetric dense 
matter. 

X^ = X^ + ^^ = y + </ , X^ = X^ + ^^ = {) + ^^ and ^ = i, + a, = (5,(t)A + Q, (3.7) 

Since the only non-vanishing background gauge field is At = A, as in Sec. |2|, we can 
put Ag = 6g(^f^A^ where (t) denotes the time component. Then the covariant derivative is 
expressed as 

D,X''^d,X'' + i[A„X''] 

=dsX'' + ds^" + + i[a., + i[as, . (3.8) 

Now we expand the action (|3.4D in terms of these fluctuations ip"". To evaluate the deter- 
minant in the action, we introduce a^s 



Grs + DrX''DsX\Gab " X'^]GacGdb) + T-^Fr^ 



(3.9) 



where 

= Grs + Gyydrydsy + T-\drAs - 9, A) , 

= Gyy{dry{ds(p^ + i[As, f^] + i[as,y]) + (9^/ + i[Ar, f^] + i[ar,y])dsy} 

+ T^^{dras - dsttr + i[Ar,as] + i[ar,As\) , 
=Gyy{idry[as,>f^] +i[ar,'/]dsy + {dr-ip^ + i[Ar , ip^] + i[ar,y]){dsP^ + i[As,p^ 

+ G^^drp'^ + i[Ar,^'>]){dsp-> + i[As,p->]) 

- iTGyyG^^{dry{dsP^ + i[As,p^]) - {drP^ + i[Ar,p^])dsy}[y,p] 



i[as,y\)} 



(3.10) 
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Then we arrive at 

V-det(a^,) =-y-det(a(0)) • det[l + (a(o))-i(a(i) + a(2))] 

= ^-det(a(0))(l + itr[(a(0))-ia«] + i(tr[(aW)-ia«])2 

- itr[((a('^))-a(^))^] + ^tr[(a(°))-a(^)] + •.•). (3.11) 
Finally, we obtain the action expanded up to second-order in the fluctuation. 
S= -^i^jd'a STv[e-y-det(a(0))(l + \ir[{a^''^)-' a^^^] + ^(tr[(a(o))-^a(^)])' 

- itr[((a(°))-^a(i))^] + ^tr[(aW)-^a(2)] + ^[X\X'^]G,,[X^,X'']G^ + ■ • •) 

(3.12) 

Now we first consider ip^ = ip which corresponds to a pseudo-scalar meson in the dual 
gauge theory. Since the fluctuation ip^ does not mix with any other fields, we could see the 
density and asymmetry dependence of the meson mass easily. In addition, the in-medium 
properties of may reveal some feature of in-medium kaon mass which is related to the 
kaon condensation in dense matter. In case wc deal with one light and one intermediate 
quarks, ip may be considered as kaons in QCD, though our model does not possess non- 
Abelian chiral symmetry. Since the mass of the strange quark is not negligible compared 
to the intrinsic QCD scale Aqcd ~ 200 MeV, the Goldstone boson nature might not be 
crucial for the mass of kaons. Now the relevant part of action is given by 

V = -/X6 y"dVS'Iv[e-y-det(a(0))(l+itr[(a(0))-ia(2)]+^[X«,X'']G;,e[X^X'^]Gd„)" , 

(3.13) 

where 

= G^^{drip + i[Ar, <p]){ds(p + i[As, ip]) 

- iTGyyG^^{dry{ds<p + i[As, <p\) - {dr<p + v])dsy}[y, ^] ■ (3.14) 

Using the explicit form of metric a^s, we evaluate the trace in the action as 

tr[(a(°))-^a(2)] = a[o>?.) = a^-^a^^ + a\^/^ + ag^aj? + a^^/^ 

1 



Gu{Gpp + Gyyiy'y} + iT-^A')^ 

X [{Gpp + Gyyiy'f}G^^iDtpf - 2iGyyG^^A'y'Dtip [y, ip] + GuG^^p'f] ■ 

(3.15) 

Here we have omitted the terms 

T-^G^^A'{^'Dtip - Dtip (3.16) 
due to the symmetrized trace. The remaining part of the action is calculated to be 

[X\X'>]G,,[X^X'']G^ = -2GyyG^^[y,pf. (3.17) 
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Now we are ready to read off the Lagrangian up to quadratic order in (p. 

^\^±{Hi±f + \^±[^±'? (3.18) 
4l2i = ^([™"HGpp + Gyy{y'f}G4.^]y^ + [FD-\G,, + Gyy{y'f}G^^]y_) 

X {ip\2^'2l + iA^{Lpi2Lp2l - <^12V'2l) + {A^f Lpi2'P2l} 

-^{[FD-^GyyG4,4,^y\^ + [FD-^GyyG^4,A'y\_Y^i2^2i - ^\2^2X - 2iA^^i2^2i) 
+]^{\FD-^GttG4,^]y^ + [FD-^GuG^^]y_)ipi2'v2i 

tHvs? f - - \ 

+ [[FGyyG^^]y^ + [FGyyG^^]y_ jipi2^21 

= ^P {ifh^Phl + iA3{(pi2'P2l - 9jl2'/i'2l) + {A3f(pi2ip2l} 



'2f„, \2 



(95129^21 - V5i2</'2i - 2iA3ipi2(p2i) + 9312 V21' + ^— </'i2¥'2i (3.19) 
Here the component fields are defined as 

if = (foT +ipiT = - \ \ = \ , (3.20) 

2 \ipi + lip2 (fO- ipsj \f21 V-J 

and the common function F as 

F = e-*-^-det(a(0)) . (3.21) 
Now we turn to gauge field fluctuations. 

a.it,p)=h ^ (3.22) 
\a2i a-y 

Like the scalar fluctuation tp^, the gauge field fluctuations decouple from the other modes. 
Therefore we can consider this mode only. Then non- vanishing a^s in ( 3.10| ) are 



4°) = Grs + GyydrydsV + T-\drAs - dsAr 

a^l^ =T-\d, + i[A,a,])=T-^Dta, 

^zt — "'tz 

Op^ = iGyyy'[az,y] + T^'^a'^ 
oip = iGyy[az,y]y' - T'^a' 



zp — "^yyl'^ZT y\y ^z 
af} = -Gyy[a,,yf . (3.23) 

Then quadratic Lagrangian becomes 

= -^ji^jd'a STv[e-y-det(a(0))(-itr[((aW)-ia(i))'] + ^tr[(a(0))-ia(2)])] . (3.24) 
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Also we can write the trace part explicitly up to overall factor 1/2 



1 



- {T^GuGpp + A'^)Gyy[y,a,f + GuKf] (3.25) 
where D = Gu{Gpp + G^^j^y'^) + (T-^A'f. Here we have also omitted the terms 

(T^DG,,)-^ [A'iDta, a', - a', Dta,) + iGuGyyy' {[y, a,]a', - a',[y, a,])] (3.26) 

under the consideration of the symmetrized trace. We can read off the Lagrangian for this 
fluctuation up to quadratic order. 

= \[F{T^DG..)-'{Gpp + Gyy{y'f}\^{d±f + \[F{T^DG..)-'Gu]y^{a±'f (3.27) 

41 = \mT^DG,.)-'{Gpp + Gyy{y'f}]y^ + [F{T^DG..)-\Gpp + Gyy{y'f}]y_) 
X {012021 + ^^3(012021 — ai2a2i) + (^3)^ai2«2i} 



-^[[F{T'DG..)-'GyyJ^yX^ + [F{T'DG..)-'GyyA'y% 

X (ai2«2i - ai2"2i - 21^3012021) 
+l([F{T''DG,,)-'GuL + [F{T''DG,,)-^GuL)au'a2i' 



+ ^y^{[F(T'DG,,)-\T^GuGpp + {A'f)Gyy]y^ 

+ [F{T^DG,,)-\T^GttGpp + {A'f)Gyy\y_yi2a2i 
= {ai2a2i + i^3(ai2a2i - 012021) + (^3)^012021} 

--|^Q (012021 - 0'i2021 - 2^^3012021) + -M0i2'02l' + ^^^^ S 012021 (3.28) 

3.2 Meson mass 

Now we calculate the pseudo-scalar and vector meson masses in asymmetric matter. We 
expect that the off-diagonal meson masses will be sensitive to the asymmetry, which is 
analogous to vr^ in isospin matter, while diagonal part is dependent only on the total 
number density like vr*^ in isospin matter. For illustration purpose, we take A = 6 and 
Mkk = 1 GeV. In Appendix |A[ we use a bit larger value for the 't Hooft coupling A = 17. 
To discuss pion-like and kaon-like mesons, we consider m+/m~ = 3 and m+/m^ = 30 in 
this section. 
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We first consider the pseudo-scalar meson which is analogous to the pion or kaon 
depending on the quark mass ratio. The equations of motion for ipi2 and (p2i are, from 
( |3l9D , 



= P(^-i2 + 2i(^3P - ysQ) ^'12 + M^'i's + EV'i2 + {-(^3)'IP + 2y3^3Q - T^ivz?^} , 

= P(^2i - 2i(^3P - ysQ) ^21 +^^21 + KV21 + {-(^3)'P + 2y3^3Q - T^y-sfS} 9931 • (3.29) 



For (p±, we have 

If we decompose ip as 



= P± V3± + M± ip'^ + R'^ip'^ . (3.30) 



<p{t,p)=e-"^'<^{p) , (3.31) 
the equations of motion can be written as 

M $;'2 + R'$;2 - [ (wi2 - ^3)'F + 2^3(^12 - ^3)Q + T^ysf S]<i>i2 = , 

M $'2^ + M'$'21 - [ i^21 + ^3)'IP - 2^3(^21 + ^3)Q + r2(y3)2S] $21 = , 

M± $'4 + M'±^>'± - w^P±$± = . (3.32) 

As expected the equations for ip± are independent of the asymmetry, which is encoded in 
^3, and depend only on the total charge Q. We solve these equations of motion numerically, 
and the results are in Fig. where the vertical dashed line is for the transition density Qc 
from a = 0, maximally asymmetric matter, to a 7^ 0. Here we choose Mkk = 1 GeV and 



A = 6 [22|. In this case, Q ~ 1.2 corresponds to the normal nuclear matter density pQ. The 



relation between the baryon number density and Q is given by 

2 • 2^/"^ 

^=81(2^^^^^^' (^-^^^ 

For illustration purpose, we show low density data only in Fig. ^(a) whose high-density 
behavior is similar to Fig. ^(b). The UV boundary condition is unambiguously fixed by 
the normalizability condition, <I> ~ 1/ p\p^^, and we choose Neumann boundary condition 
at IR. 

In Fig. ^, green lines are for bound state of two light quarks of the mass m^, and 
for (^+, bound state of two massive quarks of the mass m"*". Note that the mass of is 
independent of the density until the transition density. This can be easily understood since 
the repulsive force between the U(l) charges of two flavor branes is inactive when Q < Qc- 

For the off-diagonal part, both the total charge (or baryon number density) and the 
asymmetry will affect the meson mass. Here one can regard the asymmetry in ^3 as isospin 
chemical potential or isospin number density. In this case one can easily expect that will 
induce the mass splitting between ipi2 and ip2i since they have opposite isospin charges. 
Analogy is again vr^ with nonzero SU(2) isospin chemical potential in isospin matter: with 
increasing (m^-) goes up (down) linearly. While A3 is causing the mass splitting, 

the total charge Q would make the mass of ipi2 and ip2i decreasing as one can see from 
the behavior of the diagonal mesons, ip±. Blue (red) lines in Fig. Q are for ipi2 {<P2i), 
respectively. At low densities below Qc, A3 is dominating and the mass splitting is there. 
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(a) 



(b) 



Figure 4: The in-medium mass of the pseudo-scalar meson with A = 6 and M^k = 1 GcV. Here 
the vertical dashed line is for the transition density Qc from a — 0, maximally asymmetric matter, 
to a ^ 0. (a) For jmr = 3 (b) For jmT = 30 



Now as we increase the density Q, the asymmetry in ^3 is becoming small and so Q is 
dominating to give almost the same behavior in the mass of both off-diagonal mesons at 
high densities. 

An interesting and/or strange behavior pops up near the transition density. We focus 
on Fig. 0(b) since it shows what is going on near the transition density clearly. The mass 
of (^12 suddenly drops and that of 9921 abruptly goes up, reverting the mass hierarchy of 
9J12 and (/921- How can we understand these bizarre results? Here we try to address this 
question based on the Pauli exclusion principle, though the behavior near the transition 
in Fig. ^(b) is too drastic to be explained solely by the Pauli principle. Since the off- 
diagonal meson consists of one light (up or down) and one intermediate (strange) quarks 
for va^ jmr = 30, we may consider them as K^: (pi2 as (us) and 9921 as K~ (su). In 
Fig. §(b), the transition happens from pure up quark matter to strange matter with both 
up and strange quarks. Based on the Pauli principle between the quarks in the kaons and 
in the vacuum, we can expect that after the transition the mass of K~ would increase since 
the strange quarks are piling up in the vacuum, while the mass of drops much faster 
since the population of the up quarks reduces after the transition. This can be also seen 
in a simple model calculation, see Eq. (15) in the first arXiv version of ||2l|, where it is 
shown that the mass correction for charged kaons is proportional to i{pu — Ps)- Here pq is 
the number density for a quark flavor q. 

Though the off-diagonal meson in Fig. ^(b) is not exactly the charged kaon in nature 
partly because our D4/D6/D6 model is lack of non-Abelian chiral symmetry, we compare 
our results with the in-medium charged kaon masses obtained in a QCD effective model. 
In the context of an extended relativistic mean field model the authors of [19] studied in- 
medium kaon mass in hyperon-rich dense matter. Their results on are qualitatively in 
agreement with our results when it comes to a general tendency except near the transition 
density: m^± shows mass splitting and then both of them start to decrease at high density. 

Now we move on to the vector meson. The equations of motion are given by, from 
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(a) 



(b) 



Figure 5: The in-medium mass of the vector meson with A = 6 and M^k = 1 GcV. Here the 
vertical dashed Une is for the transition density Qc from a — 0, maximally asymmetric matter, to 
a 7^ 0. (a) For m^/m~ = 3 (b) For m~^/m~ = 30 



(IOTD and (|3j|), 

= M±^'^ + R'^A'^ - ujlf±A± , 

= M^'/2 + ^'A2 - {(^12 - AsfF + 2^3(^12 - ^3)Q + (y3)'§}-Al2 , 

= + - {{uJ2i + AsfF - 2y3(a;2i + ^3)Q + {y3f^}A2i . (3.34) 

Similar to the pseudo-scalar meson, we solve these equations numerically to obtain the 
result in Fig. ||. As it is manifest from the figure, the vector meson shows a behavior 
similar to the pseudo-scalar meson. 



4. Summary 

We studied the in-medium meson mass in a holographic asymmetric dense matter. In this 
work we considered the pseudo-scalar and vector mesons. We used the D4/D6/D6 model 
with two quark flavors to consider asymmetric dense matter. A drawback of this model 
is the absence of the non-Abelian chiral symmetry, while it possesses the vector Nf flavor 
symmetry. We studied a case with /m" ~ md/rriu to see the effect of isospin violation 
on the meson masses. Then, we took ra^ jmT ~ nis/mq, where ruq ~ [rriu + mrf)/2, to 
calculate in-medium kaon-like meson masses. In both cases we observed the mass splitting 
of charged pseudo-scalar and vector mesons at low density due to the asymmetry, while 
at high densities charged meson masses become degenerate. Near the transition density, 
we found an exotic behavior in charged meson masses. For instance, at the transition 
density the mass of Lpi2 (K^) tumbles down steeply and that of ip2i (K^) increases very 
quickly. We could partly understand this behavior based on the Pauli exclusion principle 
in a boundary quark model picture. 
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A. Mesons with A = 17 

In this Appendix we take A = 17 to see if our results are stable with a different value of 
the t' Hooft coupling constant. 




(a) (b) 

Figure 6: The in-medium mass of the pseudo-scalar meson with A = 17 and Mkk — 1 GeV. Here 
the vertical dashed line is for the transition density Qc from a — 0, maximally asymmetric matter, 
to a ^{). (a) For m'^ /m^ = 3 (b) For m'^ /m^ = 30 




Comparing with the case of A = 6, we could see that our results are not that sensitive 
to a moderate variation of the 't Hooft coupling. 
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